Oscillating axion bubbles as alternative to supermassive black holes at galactic centers 



Anatoly A. Svidzinsky 
Department of Physics, Institute for Quantum Studies, 
Texas A&M University, College Station, TX 77843-4242 
asvid@]ewel. tamu. edu 

(Dated: February 5, 2008) 

Recent observations of near-infrared and X-ray flares from Sagittarius A*, which is believed to 
be a supermassive black hole at the Galactic center, show that the source exhibits about 20-minute 
periodic variability. Here we provide arguments based on a quantitative analysis that supermassive 
objects at galactic centers may be bubbles of dark matter axions rather than black holes. An 
oscillating axion bubble can explain periodic variability of Sagittarius A* and yields the axion mass 
about 0.6 meV which fits in the open axion mass window. The bubble scenario with no other free 
parameters explains lack of supermassive "black holes" with mass M < IO'^Mq. Low -mass bubbles 
decay fast and as a result are very rare. We also found that the mass of an axion bubble can not 
exceed 1.5 x lO^M©, in agreement with the upper limit on the supermassive "black hole" mass 
obtained from observations. Our finding, if confirmed, suggests that Einstein general relativity is 
invalid for strong gravity and the gravitational fleld for the bubble effectively becomes repulsive at 
large potential. Imaging a shadow of the "black hole" at the Galactic center with VLBI in the next 
decade can distinguish between the black hole and the oscillating axion bubble scenarios. In the 
case of axion bubble, a steady shadow will not be observed. Instead, the shadow will appear and 
disappear periodically with a period of about 20 min. 



I. INTRODUCTION 

Originally introduced to explain why the strong inter- 
action, in contrast to weak interactions, does not vio- 
late CP symmetry [l| , hypothetical axions have since be- 
come one of the leading particle candidates for the cold 
dark matter in the Universe 0| . The axion appears as a 
pseudo Nambu-Goldstone boson of a spontaneously bro- 
ken Peccei-Quinn symmetry whose scale / determines 
the axion mass m, 

mw-^ = 0.62eV j (1) 

and suppresses the coupling to Standard Model parti- 
cles, (X 1//. Here ttitt = 135 MeV is the neutral pion 
mass and = 93 MeV its decay constant 0, Q . Astro- 
physical and cosmological arguments constrain the axion 
mass m to be in the range of 10^^ — 3 x 10"'^ eV 0. 
Axions in this mass range could provide much or all of 
the cold dark matter in the Universe. Properties of stars 
impose the upper limit on the axion mass via constraints 
on axion interaction with photons, leptons and nucleons. 
However, such interactions are model-dependent. Obser- 
vations of the cosmological large-scale structure constrain 
the axion-pion coupling which yield weaker upper mass 
limit m < 1.05 eV [1]. 

Interaction of axions with QCD instantons generates 
the axion mass and periodic interaction potential Q 

yM=m2/2[l-cos(^//)], (2) 

where (/? is a real scalar axion field. 

Here we argue that oscillating axion bubbles, rather 
then supermassive black holes, could be present at galac- 
tic centers. Recent observations of near-infrared and X- 
ray flares from Sagittarius A*, which is believed to be 



a 3.6 X 1O^M0 black hole at the Galactic center, show 
that the source exhibits about 20-minute periodic vari- 
ability [1, 0, [1] ■ An oscillating axion bubble can explain 
periodic variability of Sagittarius A* and yields the axion 
mass about 0.6 meV. Fig. [T] explains the mechanism of 
hare variability. 

Moreover, the bubble scenario explains observed lack 
of supermassive "black holes" with mass M < IO^Mq. 
As we discuss later in this paper the bubble life-time 
t (X M^/^, it becomes less then the age of the Universe for 
M < 5x lO^M©. The bubble at our Galactic center would 
decay within about 5 x 10^ yrs. If, however, M < lO^M© 
the decay time becomes very short, t < 10^ yrs, and as a 
result such objects are very rare. 

Finally, the axion bubbles with no free parameters (if 
we fix m = 0.6 meV based on Sagittarius A* flare vari- 
ability) explain the upper limit (1.5 x lO^M©) on the 
supermassive "black hole" mass found in recent analysis 
of the measured mass distribution 0. 

In recent years, the evidence for the existence of 
an ultra-compact concentration of dark mass associated 
with the radio source Sagittarius A* in the Galactic Cen- 
ter has become very strong. However, an unambigu- 
ous proof that this object is a black hole is still lack- 
ing. A defining characteristic of a black hole is the event 
horizon. To a distant observer, the event horizon casts 
a relatively large "shadow" with an apparent diameter 
of about 10 gravitational radii due to bending of light. 
The predicted size (^--^30 micro-arcseconds) of this shadow 
for Sagittarius A* approaches the resolution of current 
radio-interferometers. Hence, there exists a realistic ex- 
pectation of imaging the shadow of a black hole with 
very long-baseline interferometry within the next decade 
[Tol ITTl Ha . [l3| . Such imaging will allow us to distinguish 
between the black hole and the oscillating axion bubble 
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axion bubble 




is the dimensionless pot ential and the couphng parameter 
respectively, mpi = yJhc/G — 1.2 x 10^^ GeV is the 
Planck mass. 

The interaction potential V has degenerate minima 
V — Q a.t (p — 27rn/a, where n is an integer number. As 
a consequence, equation for the axion field </? has bubble- 
like solutions. The bubble surface is an interface between 
two degenerate vacuum states with ip = 27rn/a (r < R) 
and If = Q {r > R). In this paper we consider spherical 
bubbles with surface width much smaller then its radius 
R and n ~ 1. Energy density of the axion field is nonzero 
only at the bubble surface. Energy of a static thin-wall 
bubble is = A-kcjR^, where a is the surface energy per 
unit area (which equals to the surface tension for the do- 
main wall we consider Il4|) determined by an integral 
over one potential period jisl. [l6| 



(6) 



FIG. 1: Mechanism of flare variability. Axion bubble radius 
oscillates between -Rmin and i?max which yields periodic vari- 
ation of the gravitational redshift z{t) of the bubble interior. 
The dark matter bubble itself does not radiate electromag- 
netic waves; radiation is produced by flaring clouds of baryon 
matter trapped in the bubble interior. Redshift z{t) reduces 
radiation power of the cloud by a factor of 1/(1 -I- As 
a result, intensity of the flaring cloud radiation is modulated 
with the period of bubble oscillation. Modulation amplitude 
depends on the distance between the cloud and the bubble 
center. 



scenario which we propose in this paper. If the axion bub- 
ble, rather then a black hole, is present at the Galactic 
center, the steady shadow will not be observed. Instead, 
the shadow will appear and disappear periodically with 
a period of about 20 min. 



II. AXION BUBBLES 

We introduce dimensionless coordinates and define the 
unit of distance, time and ip as 

h h 1 



ro = 



mc 



(3) 



where c is the speed of light and G is the gravitational 
constant. For the moment we omit gravity. Further we 
use natural units for which h = c = 1. Energy of the 
axion field in units of ni^^/m is given by 



E 



dr 



where 



V{cp) 



:[1 - cos(a(^)]. 



a 



(4) 



(5) 



Surface tension a depends only on the axion interaction 
strength. The later, however, slightly depends on the 
axion model which can change Eq. ([6]) by a factor of the 
order of one [l3l- 

Under the infiuence of surface tension an initially static 
bubble collapses to its center which yields reduction of 
the surface energy. However, if we include gravity this 
gives an additional energy contribution. Such a contribu- 
tion could substantially alter the bubble evolution and, 
in particular, prevent collapse as we discuss in the next 
section. 

One should mention that decay of axions into photons 
is suppressed in the bubble. Such a decay is not allowed 
by energy conservation. The bubble surface is approxi- 
mately a one dimensional kink. In the kink's reference 
frame the kink is static and the distribution of the ax- 
ion field is obtained by minimization of the energy func- 
tional (HI) subject to the boundary conditions that far 
from the kink we have fixed vacuum states ((/? = from 
one side and ip — 2T:/a from the opposite side). The op- 
timized field distribution determines the kink's energy. 
Any small change in ip would increase the total energy 
of the axion field. Axions in the kink cannot decay into 
photons because annihilation of the axion would change 
the distribution of the axion field in the kink and, hence, 
increase the kink's energy. In such a process both the ax- 
ion field and photon acquire energy which violates energy 
conservation. 



III. THE ALTERNATIVE THEORY OF 
GRAVITY VS EINSTEIN GENERAL 
RELATIVITY 

So far Einstein general relativity has successfully 
passed all available tests. However such tests have in- 
spected the theory only at weak gravitational field p^ . 
One should note that observations of binary pulsars 
yet have not provided a test of general relativity at 
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strong gravity. Rather, such observations tested Einstein 
equations in the post-Newtonian approximation and the 
strong equivalence principle [H]. 

Are Einstein equations valid for strong gravity? The 
answer to this question remains unknown and only ap- 
propriate observational tests can shed light on it. It is 
well known that in Einstein general relativity the gravita- 
tional field disobeys the principle of superposition. This 
is the consequence of the postulate that the space-time 
metric is determined by the Einstein equations 

1 , , 

Rik - -jpikR = SvrTifc, (7) 

with the particular choice of the matter energy- 
momentum tensor T^^ proposed by Einstein. However 
the Einstein theory can be modified by modifying the 
energy-momentum tensor. This yields a possibility to 
satisfy the superposition principle by a proper choice of 

In Appendix A we derive a space-time metric produced 
by a static mass distribution based only on the superpo- 
sition principle. The answer is given by the Yilmaz ex- 
ponential metric. Then we show that the Einstein equa- 
tions yield the exponential metric if Tifc is taken as an 
"electrostatic" energy-momentum tensor. 

In the weak field limit the exponential metric is equiva- 
lent to those obtained in Einstein theory and, hence, the 
exponential metric agrees with the four classic tests of 
general relativity. However in the opposite limit of strong 
gravity the exponential metric is dramatically different. 
Since the superposition principle is satisfied the expo- 
nential metric predicts no black holes, but rather stable 
compact objects with no event horizon and very large, 
but finite, gravitational redshift ("dark red holes"). This 
suggests that gravitational field for those objects effec- 
tively becomes repulsive at large gravitational potential. 

Here we analyze properties of compact objects at galac- 
tic centers and show that they are in favour of the expo- 
nential metric. Our conclusion is based on a quantitative 
analysis which is independent of the particular choice of 
the time-dependent theory of gravity [li|. This is pos- 
sible because the main part of the bubble dynamics we 
use for the quantitative comparison occurs in the well- 
tested limit of Newtonian gravity. Only a small part of 
the trajectory near the lower radius turning point (where 
gravity effectively becomes repulsive) is beyond Newto- 
nian description. As a result, e.g., the period of bubble 
oscillation can be accurately obtained using Newtonian 
gravity, independent of which theory of gravity yields the 
repulsive force at small radius. 



A. Bubbles in exponential metric 

In Appendix A we derive a metric for a static mass 
distribution assuming that the gravitational field obeys 
the principle of superposition for any field strength. The 



answer is given by an exponential isotropic line element 
of the class proposed by Yilmaz [13, HH 

ds^ = -^e^'^dt^ + e-^'l'{dx^ + dy^ + dz^), (8) 

where for a spherically symmetric thin- wall bubble of ra- 
dius R (units are c — G ~ 1) 

and M is the bubble mass. The exponential metric does 
not have the singularity of the Schwarzschild solution at 
finite radius, and therefore replaces the concept of black 
holes with that of "dark red holes" . 

In the reference frame of a distant observer the energy 
of a static bubble is given by [13] 

U{R) = AncrR^ exp , (10) 

where a is the intrinsic surface energy density (as it would 
appear to an observer located at the bubble surface) 
given by Eq. ([6]), M is the dimensionless bubble mass 
in units of rn^^/m and R is the dimensionless bubble ra- 
dius in units of tq. If the radius of a spherical bubble 
changes with time then the total energy is 

E = U{R) + Ek, (11) 

where Ek > is the kinetic energy. The total energy E 
is a constant of motion, then based on the equivalence 
principle we obtain M ^ E =const. Evolution of the 
bubble radius is similar to a one dimensional motion of 
a particle in the effective potential U{R). We plot U{R) 
in Fig. [21 The effective potential has a shape of a well 
and depends on the total energy (mass). At i? ^ M one 
can omit gravity and U{R) ~ AnaR"^ is just a surface 
energy (tension) which tends to contract the bubble. At 
R <C M gravity produces large repulsive effective poten- 
tial which forces the bubble to expand. As a result the 
bubble radius R{t) oscillates between two turning points 
determined by U{R) = M. 

Fig. [3] shows numerical solution of the equation for 
turning points M = AiraR'^ exp {M/ R). For AiraM < 
4/e^ = 0.541 equation has two solutions; radius of such 
bubbles oscillates with time between turning points Ri 
and Rf). If M — A/max — l/{e'^na) the bubble is static 
with R = M/2. Such static bubble possesses maximum 
possible mass. Bubbles with M > iWmax do not exist. 
As we show below, for the Galactic center bubble AI ^ 

In a general case to describe R{t) quantitatively we 
need to use dynamic equations. These equations depend 
on a particular choice of the time-dependent theory of 
gravity. In this paper, however, we do not need them 
which makes our results quite general. The point is that 
if the Galactic center object is an axion bubble the main 
part of its periodic oscillation occurs in the limit R{t) 3> 
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FIG. 2: Effective potential for axion bubble motion U{R) 
(solid line) in Yilmaz exponential metric. Bubble radius R{t) 
oscillates between turning points determined by the equation 
U{R) = E. We expressed R in Astronomical Units (AU) and 
U (R) in solar masses, m is the axion mass. 




FIG. 3: Mass-radius relation for a bubble in Yilmaz expo- 
nential metric. At a given mass A/ the bubble radius R{t) 
oscillates between turning points Ri and -Ro- If M = Mmax 
then Ri = Rq and the bubble is static. 



M. In this region we can omit gravity and use the well- 
tested special theory of relativity that yields the following 
mass-radius equation for a relativistic bubble 



M 



(12) 



^1 - [dR/dty 

which has a simple solution 
where Rq = TTcr is the maximum bubble radius and 



(13) 



cn(a;, k) is Jacobian elliptic cosine. For small R{t) the 
solution (|13p is not applicable. In this region the bubble 
shrinking slows down and after reaching the inner turning 
point i?i the bubble starts to expand. Assuming that Eq. 
p3|) is accurate for the main part of the motion we obtain 
that the period of bubble oscillation is T ~ 2.622i?o- 
Taking into account that M = AnaB^ = 8Rq /a^ , we get 



in dimension units T 
Eq. (P) we find 



0.262^/ M/mh/fc^. Then using 



r~ 0.523 



M 



lOSMm 10-3eV 



X 15.27 min. 



(14) 
8] then Eq. 
IQio GeV). 



If M = 3.6 X lO'^M© and T = 22.2 min 
(fn)l yields the axion mass m ~ 0.6 meV (/ 
One should mention, however, that due to time dilation 
the period of flare variability depends on the distance 
between the flare source and the bubble center (see Fig. 
[IJ and, thus, could differ from Eq. by a factor of 

the order of one. This yields an inaccuracy in the axion 
mass determination in the same factor. 

Next we discuss the bubble life time. The bubble de- 
cay occurs by means of axion emission. Due to spheri- 
cal symmetry there is no radiation of gravitation waves. 
Bubble surface, the interface between different vacuum 
states, is a soliton (or a kink) that is studied in many 
areas of nonlinear physics. One dimension solitons, con- 
trary to 2D or 3D, are stable and preserve their shape 
under reflection from a boundary. Because a thin-wall 
bubble surface can be treated as a ID soliton this insures 
its very long life time. However, due to finite bubble ra- 
dius the ID treatment is only approximate. Deviation of 
the problem from ID leads to slow decay of the soliton 
by emission of particles (axions). 

We estimate the bubble decay rate as the time of en- 
ergy loss by the bubble with the radius R{t) oscillating 
between the outer Rq and inner i?i turning points. En- 
ergy loss by the bubble surface becomes substantial only 
when R{t) < Rq^^ (see Appendix B below). In our case 

Ri 3> Rq and, therefore, the region of intensive en- 
ergy dissipation is not accessible. As a result, the energy 
emission is negligible yielding long-lived bubbles. In Ap- 
pendix B we estimate the bubble life-time. The answer 
is given by Eq. (|B11[) which in dimension units reads 



t ^ 



Rq ( Ri 

Rq 



Ri 
ro 



(15) 



If TO = 0.6 meV then the bubble surface has thickness 
^0 — h/mc = 0.3 mm. For a bubble with mass M = 3.6 x 
1O^M0 the maximum radius is Rq = A83Rq, while the 
gravitational radius i?g = IG.IRq. One can find the inner 
turning point Ri from the equation Rq — R\ exp(A//i?i) 
which yields Ri — 1.\Rq. Using Eq. we then obtain 
the bubble life time i ~ 5 x 10^ yrs. This can explain 
observed lack of supermassive "black holes" with M < 
1O^M0. Axion bubbles with such masses decay fast with 
life time t cx M^/^. 
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For bubbles with M <C Mmax one can obtain Ri in 
terms of the bubble mass M analytically 



M 



M 




(16) 



Substitute this into Eq. yields the following expres- 
sion for the bubble life-time 

^ V8c2mi2ln6[a;ln(a;)] ~ 



M 



(i06Mq 



9/2 



4.73x10" 
X — — —T years, (17) 



10-3eV ln^[a;ln(x)] 



where 



X = 



'pi 



277771 ,r/7rVM 



140 



m lO^Af© 
10-3eV M 



Eq. dni) shows that the bubble life-time t cx M^/^ ^j^j 
it becomes less then the age of the Universe for M < 
5 X IO'^Mq. For M < IO^Mq the decay time becomes 
very short, t < 10^ yrs, this is why we do not observe 
supermassive "black holes" with such masses. 

Recent analysis of the mass distribution for the com- 
pact objects at galactic centers shows existence of an up- 
per limit for the supermassive "black hole" mass 



Mn 



1.2tli X IO^Mq. 



(18) 



Next we calculate the maximum possible mass of an axion 
bubble. In dimension units it is given by 



M„ 



(19) 



For m = 0.6 meV Eq. HH) yields M^ax = 1-5 x IO^Mq. 
This value agrees with the upper limit on the supermas- 
sive "black hole" mass measured for galactic nuclei. 
Radius of the static bubble with Mmax is i? = 1673i?0. 



B. Bubble in Einstein general relativity 

In Einstein general relativity for a spherically symmet- 
ric bubble the metric can be written in the form 



-hdt^ + gdr^ + r^dVP 



(20) 



where the radial metric, and h, the lapse, are functions 
of t and r with r being the circumferential radius. For a 
static thin- wall bubble of radius R the Einstein equations 
yield 



h{r) 



1 - 2M/r, r ^ i? 
1 - 2M/R, r <R, 



(21) 




E=const 
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FIG. 4: Effective potential for axion bubble motion in Ein- 
stein general relativity. "Small" bubbles shrink to black holes, 
while large bubbles expand infinitely. We expressed 7? in As- 
tronomical Units (AU) and U{R) in solar masses. 



q(r) = <! 1-2MA' > ^ 
' 1 1, r<R. 



(22) 



Note that g{r) undergoes a jump at the bubble surface, 
while exponential metric ([8]), ([9]) is continuos. Energy of 
the thin- wall bubble with radius R{t) is given by [U, [2^ 



E 



inaR'^ 



y/l - {dR/dry 



(23) 



where r is the interior coordinate time (dr^ = hdt^). The 
corresponding effective potential 



U{R) ^ AnaR^ 



iir^a^R^ 



(24) 



is pictured in Fig. [4l 

In Einstein theory "small" bubbles shrink toward the 
gravitation radius Rg = 2M and at t 3> Rg/c behave as 
black holes, while large bubbles expand infinitely. This is 
dramatically different from bubble evolution in the expo- 
nential metric which we discussed in the previous section. 
At the same time, the exponential and the isotropic form 
of the Schwarzschild metric (jA2[) are the same to second 
order in the gravitational potential in the temporal 
part and to first order in the spatial part (see Appendix 
A). This is sufficient to insure that they both give iden- 
tical results in the four classic weak-field tests of general 
relativity. 



IV. DISCUSSION 

If axion bubbles, rather then supermassive black holes, 
are located at galactic centers then what is the mecha- 
nism of their nucleation? Dark matter axions, if they 
exist, form halos around galaxies. The halo of axions is 
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in a quantum degenerate non-equilibrium regime. Evo- 
lution of the axion halo is governed by the self-gravity 
and axion interaction V{ip). The interaction V{ip) be- 
comes important only for dense axion clumps. Dynam- 
ics of a dilute galactic halo is determined by self-gravity. 
Seidel and Suen have studied evolution of a massive, self- 
gravitating real scalar field in Newtonian limit (omitting 
self- interaction V{(p)) [251. They have shown that in- 
dependent of the initial conditions a scalar field config- 
uration collapses to form a compact object by ejecting 
part of the scalar field, carrying out the excess kinetic 
energy. The cooling occurs due to nonlinear effects of 
the self-gravitation of the field. Characteristic cooling 
time is a free falling time to the center due to self grav- 
ity t ~ 2i?j^{f^/V2GMhaio, where i?haio and Mhaio is an 
initial radius and mass of the axion halo in a galaxy. For 
^haio = 60 kpc and Mhaio = lO^^M© we obtain the char- 
acteristic cooling time t ^ 10^ yrs. 

Thus, within about 10^ yrs the gravitational cooling 
mechanism yields formation of compact axion clumps in 
a galactic halo. Evolution of such clumps is then gov- 
erned by self-interaction V{(p) which leads to bubble for- 
mation. This is shown by three-dimensional numerical 
simulation of the evolution of inhomogeneities in the ax- 
ion field due to the self-interaction V{(f) [l^. Such a sim- 
ulation (which omits gravity) has indeed demonstrated 
formation of bubble-like structures (see Fig. 5a in Ref. 
[HI ) . The mass of the nucleated bubbles is much smaller 
then the mass of the halo they are born in. However for 
typical halos many bubbles are born with masses much 
greater then IO^Mq and hence they are long-lived ob- 
jects. 

In Einstein general relativity axion bubbles under the 
influence of surface tension collapse fast into black holes. 
However, the Einstein theory yet to be tested in the limit 
of strong gravitational field. There is a possibility that 
at strong field the gravity is not described by Einstein 
general relativity, and rather by an alternative theory 
which also passes all available tests. In this paper we 
consider axion bubbles in a very general approach avoid- 
ing a particular choice of the alternative time-dependent 
theory of gravity. Our results are valid for any metric 
theory which in based on the principle of superposition. 
This principle yields the exponential metric in the static 
limit, as shown in Appendix A. One should mention that 
if the space-time geometry is described by the exponen- 
tial metric then compact supermassive objects at galactic 
centers can not be made of baryonic matter. Maximum 
mass of a compact (neutron star like) baryonic object in 
such a metric can not exceed about 12M0 [23|. Hence, 
dark matter of non baryonic origin is the only alternative 
for their composition. 

We found that in the exponential metric the axion bub- 
bles with M > IO^Mq are very long lived. Instead of 
collapsing into a black hole the bubble radius oscillates 
between two turning points determined by the net mass. 
Such oscillating bubbles, rather then supermassive black 
holes, could be present at galactic centers. Our result can 



account for periodic variability observed in near-infrared 
and X-ray flares from Sagittarius A* 0, Q and yields 
the axion mass about 0.6 meV. 

Moreover, the bubble scenario with no free parameters 
(if we fix m = 0.6 meV based on Sagittarius A* flare 
variability) explains lack of supermassive "black holes" 
with M < IO^Mq. We find that if M < IO^Mq the 
bubble life time becomes very short, t < 10^ yrs, and 
as a result such objects are very rare. We also found 
that for the exponential metric the bubble mass can not 
exceed Mmax = 1-5 x IO^Mq. This, again with no free 
parameters, explains the upper limit on the supermassive 
"black hole" mass measured for galactic centers [H . 

For axion with mass to = 0.6 meV the axion-photon 
coupling constant is ga-y ~ 10~^'^GeV~^ [1^. Recently it 
was argued that the solar corona A"— ray emission can be 
explained by solar axions of the Kaluza-Klein type (that 
is by axions propagating into extra dimensions) which 
are gravitationally trapped by the Sun and decay near 
the solar surface [13, |3l|. The estimated value of ga-y 
from the analysis of solar corona A— rays is similar to 
our finding; this is an interesting coincidence. 

Observation of the Galactic center with very long- 
baseline interferometry within the next few years will 
be capable to test theories of gravitation in the strong 
field limit. Such an observation will allow us to distin- 
guish between the black hole (predicted by Einstein gen- 
eral relativity) and the oscillating axion bubble scenario 
which we propose in this paper. If future observations 
indeed discover periodic appearance of the shadow from 
the Galactic center object this will also be a strong evi- 
dence for the axion nature of dark matter and will lead 
to an accurate measurement of the axion mass. 

I am very grateful to E. Sezgin and N. Suntzeff for 
useful remarks. 



APPENDIX A: DERIVATION OF THE STATIC 
EXPONENTIAL METRIC FROM THE 
PRINCIPLE OF SUPERPOSITION 

Let us consider a point mass M located at r = 0. Static 
gravitational field produced by the point mass possesses 
spherical symmetry. Without loss of generality one can 
look for the metric in an isotropic form 

ds^ = -h{r)dt^ + g{r){dr^ +r^d9^ + sin^ edifi'^). (AI) 

Einstein equations yield well known Schwarzschild solu- 
tion [Hi 

Here we derive h{r) and g{r) in a different way. First 
we note that a static gravitational field of any strength 
has a potential [11] 

<^(r) = In v/%) (A3) 
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and the gravitational force acting on a test rest particle 
with mass m is 



-mV0 fa 



In Minkowski space-time the potential (f> satisfies the 
Poisson equation A0 — 47rAf(5(r). Writing the Laplacian 
and the delta-function in curvilinear coordinates with 
metric gik the Poisson equation yields 



1 



d 



\/-\9tk\i 



ik 



An 



\/-\9ik\ 



MS{r) 



(A5) 

where \gik\ is determinant of the metric tensor and g*'^ is 
the tensor reciprocal to gik, that is giug^^ = 5\- For the 
metric (|Al|) Eq. (|A5jl reduces to 



(A6) 



Eq. (|A6p (with (j) from Eq. (|A3[) ) describes a relation 
between the functions h and f/ which the metric must 
satisfy. We note that Eq. (|A6[) is consistent with Ein- 
stein general relativity because the Schwarzschild solu- 
tion Jll]) obeys Eq. 

To find the functions h and g we need an additional 
constraint. Here we postulate that the force of gravity 
must obey the principle of superposition at any strength 
of the gravitational field. This is the only difference from 
Einstein general relativity we introduce. Such way of 
thoughts makes an appealing connection with the quan- 
tum theory. 

The principle of superposition is formulated in coor- 
dinate systems in which we measure space coordinates 
by ideal rods unaffected by gravity (e.g. of atomic con- 
stitution). This assures that the coordinate system is 
independent of the position and the value of masses. Let 
us consider masses Mi and M2 located at coordinates ri 
and V2 in the mentioned above coordinate system. The 
superposition principle means that the gravitational force 
on any test particle due to masses Mi and M2 equals to 
the vector sum of the force due to the mass Mi located 
at ri if there is no mass M2 and the force due to the mass 
M2 located at Y2 if there is no mass Mi . 

Please note that we formulate the superposition princi- 
ple for the covariant force vector fa as in Eq. (jA4p . This 
makes a proper connection with the Newtonian limit in 
which the measurable gravitational force is a covariant 
vector (under coordinate transformation it transforms 
like derivatives of a scalar). The covariant force fa on 
a test particle is defined as [28] — mcPxa/ds"^ , where 
Xa is the coordinate of the test particle with mass m and 
ds is the interval. If at the moment when the force is 
measured the particle has zero velocity (that is when Eq. 
(|A4[) applies) then ds = dr, where t is the proper time 
and, hence, fa — rad'^Xa/ dr'^ . This equation shows that 
the definition of the covariant force is unaffected by grav- 
ity (if the coordinates Xa are unaffected) and therefore 



fa is the relevant vector to formulate the superposition 
principle. On the other hand, the contravariant vector 
/" — 9"^ fi3 contains metric in its definition and hence 
cannot be used in the superposition principle. 

It follows from Eqs. (K^i and that the force of 
gravity satisfies the principle of superposition if and only 
if 



h{r)g{r) = 1, 
which yields for arbitrary field strength 



d 



,dct>{r) 



47rA/(5(r). 



(A7) 



(AS) 



dr \ dr 
Hence 4>{r) — ~-M/r and 

h{r) = cxp(-2Af/r), g{r) = exp(2Af/r). (A9) 



Eq. (jA9|) is known as Yilmaz exponential metric 20|. 

For small A4 /r both the Schwarzschild (|A2|) and the 
exponential (jA9p metrics yield the same expansion 

2M 2Af2 2M 
h{r) = 1 - — -1- — ^ + . . . , g{r) = 1 + _ + . . . 

(AlO) 

and hence both metrics pass the four classic weak-field 
tests of general relativity. Accuracy of current tests is yet 
far from ability to check the next terms in the expansion 
(jAlOp where the two metrics start to deviate from each 
other [18j]. 

The principle of superposition allows us to find the 
metric in the case of N point masses A/i, A/jv lo- 
cated at ri, ... tn. To apply the principle of superposi- 
tion we must choose the coordinate system in which we 
measure space coordinates by ideal rods unaffected by 
gravity. Since the speed of light measured by such rods 
and clocks (e.g. of atomic constitution) is independent 
of direction in a gravitational field the metric in such 
coordinates is isotropic, that is 

ds^ = -h{v)dt^ + g{v){dx^ + dy^ + dz^). (All) 

For the metric (|XTT|) Eq. (|^ yields 



dx 



dx 



dy 



dy 



dz 



dz 



= 47r[A/i<5(r - n) + . . . + Mn5{v - tn)]. (A12) 

The principle of superposition for the potential <j) (and 
hence for the force of gravity f ) is satisfied provided hg = 
1. Then Eqs. (jAnj) and jAS]) give the following answer 
for A^— body space-time geometry in isotropic Cartesian 
coordinates 



ds' 



ca<l>de ^ 6-2"^ (dx^ + dy^ + dz'), 



where </> is the iV— body potential 

Ml Mm 



(r) 



(A13) 



(A14) 



|r-ri| |r-rAr| 

Thus, finding the space-time geometry reduces to a sim- 
ple "electrostatic" problem. 

Next we discuss how to obtain the exponential metric 
from Einstein equations. 
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1. Derivation of exponential metric from Einstein 
equations 

Einstein equations read 

R^k - ^9^kR = 87rT,fc, (A15) 

where Rik is the Ricci tensor, R is the scalar space curva- 
ture and Tik is the energy-momentum tensor of matter. 
Let us consider N fixed point masses Mi, M]\f located 
at ri, ... rjv- For such a system the energy- momentum 
tensor in Einstein general relativity is 

N 

^^MjS{r-rj), ah other i;'^ = 0. (Af6) 

Solution of Eq. (jA15|) with this Ti^ yields black holes 
and no superposition principle. 

To obtain exponential metric we must use another 
energy-momentum tensor. Let us write Tik by analogy 
with electrostatic. For N fixed point electric charges qi, 
qn located at ri, ... rjy the energy-momentum tensor 
is given by [28j (in curvilinear coordinates) 

I 



where (V^)^ = EaE" = g°'f^{d(j)/dx°'){d(j)/dx'^) and (j) is 
the electric potential satisfying the Poisson equation 



^ -4vr^g,5(r-r,). (Af8) 

\ / 1=1 



We assume that for the system of N fixed point masses 
Ml, Mn the energy- momentum tensor Tik is given 
by Eqs. (|A17|) and (jAlSp with the change Mj. 
Substituting this tensor into Einstein equations (|Af 5[) 
we obtain the solution for the metric given by formulas 
(IIl3)l and (Km . 

One should mention that we can find the proper 
energy-momentum tensor (|Af 7p simply by plugging the 
exponential metric (jA13|) into the left hand side of Ein- 
stein equations (jAlsp . Then the "electrostatic" energy- 
momentum tensor (jA17p is obtained automatically. 

The result discussed here is valid for a static gravi- 
tational field. How to generalize it for time-dependent 
fields is beyond the scope of the present paper. 



APPENDIX B: ENERGY EMISSION FROM A SHRINKING BUBBLE 



Here we calculate energy loss by a shrinking spherically symmetric bubble caused by emission of scalar particles 
(axions). For an order of magnitude estimate one can omit the effect of gravity. Then the evolution of the scalar field 
ip{t, r) is described by sine-Gordon equation 

>f-Lp" + - sin{aip) = 2ip'/r, (Bf ) 

a 



where r is the radial coordinate. Without right-hand side, Eq. (IBf [) has an exact, so-called kink, solution 

4 f \ (r-vt-Ro)]] 

ipQ = — arctan < exp ± ■j^=^= — > , (B2) 

a I [ VI - J J 

where i?o 3> f is the initial bubble radius. The solution describes a kink (space region where ip changes from 2?:/ a to 
0) propagating with constant velocity v; the kink's size is I ^ Vf — w^. 

If I <^ R{t), where R{t) is the bubble radius, r.h.s. of (|Bip may be treated as a small perturbation. Eq. (jBip 
possesses approximate solution in the form of the kink (jB2p with parameters slowly changing in time under the action 
of the perturbation. In particular, the kink shrinks due to its surface tension so that the bubble radius and the 
velocity evolve as ^] 



R{t) = Rocn{V2t/Ro, 1/V2), v{t) = ^ 1 - R^{t)lRt (B3) 

where cn stands for the elliptic cosine with the modulus I / \/2- Such a process is accompanied by emission of scalar 
particles which yields the energy loss. We estimate the energy loss following the original work of Malomed [H, [s^. 
In terms of the inverse scattering technique, the spectral density of the emitted energy Eg^(t, q) is 

^-E^^A_\B{t,q)\\ (B4) 
aq na 
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where q is the radiation wavenumber and the perturbation-induced evolution equation for the complex amplitude 
B{t,q) is given by [H, [11] 



dB_ 
lit 



2(A2+72) J_ 



dr[y-Y- 2iA7tanh 



r ^ vt 



cxp 1 + q'^t — iqr^ dr^o, 



(B5) 



where A = ^/l + q^ — q and 7 = (1 + v)/2\/l — v"^. Calculating the integral in (IB5[) yields 

dB _ m [A2(1 - v){\ - VT+v) - v/2\ cxp {iy/l + qH - iqvtj 
Itt ^ (l + ?;)/4 + A2(l-t;) cosh [7rgVT~^/2] 

If V slowly varies with time one can take v « const in Eq. (jB6[) . then after integration we obtain 

TT [A2(l -v){l- yrr^) - v/2\ exp {i^l + qH - iqvtj - 1 



(B6) 



Therefore 



Bit,q) 



[(1 + v)lA + A2(l - v)] (v/TTg2 - qv) cosh [vrgVl - 



(B7) 



dq a2[(l + t-)/4 + A2(l-i;)]2(^ + g2 — (^77)2 cosh^ [tt^VI — t;2/2] 



(B8) 



Integration of (|B8p over gives the emitted energy as a function of time Ef,{t) = j^^o dq{dEf,/dq). In Eq. (jB8|) sine 
is a fast oscillating function, so we substitute sin^(a::) -^1/2. The radiation power increases when the kink's velocity 
V approaches the speed of light c — 1. Assuming 1 — u <C 1, integration of Eq. (jBSp yields 



Ee{t) 



2.51 



7.10 



i2(l-w(t))3/2 Q,2 



(B9) 



The emitted energy becomes comparable with the initial bubble energy Eq ~ AiraRl = SRq/o^ when the bubble 

radius reaches the value i?* « Rq^^. This value agrees with those obtained in [34| . 

If the bubble shrinks from the outer turning point Rq to the inner turning point i?i 3> i?* the radiated energy is 



E, 



7.10 /i?, 



a2 V-Ri 



(BIO) 



To emit all its energy the bubble must oscillate between Ri and i?o about Eo/Ee cycles. As a result, the bubble 
life-time is 



f 7? ^0 ^ 



(Bll) 
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